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One idea to explain the mysterious dark energy which appears to pervade the Universe is that 
it is due to a network of domain walls which has frozen into some kind of static configuration, 
akin to a soap film. Such models predict an equation of state with w = P/p = — 2/3 and can 
be represented in cosmological perturbation theory by an elastic medium with rigidity and 
a relativistic sound speed. An important question is whether such a network can be created 
from random initial conditions. We consider various models which allow the formation of 
domain walls, and present results from an extensive set of numerical investigations. The idea 
is to give a mechanism which prevents the natural propensity of domain walls to collapse and 
lose energy, almost to the point where a domain wall network freezes in. We show that when 
domain walls couple to a field with a conserved charge, there is a parameter range for which 
charge condenses onto the domain walls, providing a resistive force to the otherwise natural 
collapse of the walls. 



1 Introduction and motivation 

Modern cosmology has found itself a problem. If observational data are understood within the 
framework of the FRW model in General Relativity, one must invent a "dark energy" component 
within the content of the universe, to explain acceleration (in the absence of a cosmological 
constant). The alternative is to move away from the FRW model, or to modify the gravitational 
theory. 

In this paper we provide fresh theoretical evidence for the viability of the elastic dark energy 
model - a model which uses frozen domain wall networks to accelerate the universe^. The 
equation of state for a gas of domain walls moving at velocity v is Wdw = v 2 — 2/3, hence our 
requirement for a frustrated "frozen" network, for which = —2/3. We present a mechanism 
which freezes domain wall networks. If symmetry currents interact with domain walls, a pseudo- 
stable glass- like network forms. 

We will begin by giving a quick review of features and results of the proto-typical domain 
wall forming model. We will then discuss two models, which have the feature that a field with 
discrete minima interacts with a second field with continuous U{1) symmetry. 



2 Z.2 domain walls 



If a field theory has a potential, and if the potential has multiple disconnected vacua, domain 
walls form. Domain walls are the interpolations between spatially adjacent vacuum states^. 
The simplest example of such a field theory is described by the Lagrangian 

C= l -d^4>-\{^-r, 2 f . (1) 

The vacuum in this theory - the configurations of the scalar field ^{x^ 1 ) G M that minimise 
the potential - is just the set V = {+rj, — rf\. If we imagine two chunks of space, one with 
4> = +7] and another with <j) = —rj, then the field must continuously interpolate at the boundary 
between the "domains". There is a static solution to the equations of motion which captures 
such an interpolation, 4>(x) = tanh(x/A), where A := y / 2/\r] 2 , the so-called kink solution. 
The symmetry properties of this model are rather interesting. Consider the full model ([!]), 



and transforming <f> —> t ne model is totally invariant under this transformation. However, 
if one expands about each of the minima, the symmetry is lost. This is a simple example of 
spontaneous symmetry breaking. 

By a rather simple dimensional argument, as well as substantial numerical investigations^^ 
the number of domain walls formed by the standard Kibble scenario drops, n oc t , This result 
is also observed in many other models (so long as discrete vacua are present) - even with 3 or 
more vacua (with junctions). The standard result, therefore, is that domain walls collapse under 
their own tension. 

Numerical simulation Throughout the rest of the paper we will refer to simulations that 
we have performed. We evolve the equations of motion in (2 + l)-dimensions, using a leapfrog 
evolver, imposing periodic boundary conditions. We use random initial conditions: each grid- 
point is randomly assigned to be in one of the vacuum states of the theory. We also apply 
damping for a small amount of the total simulation time, to ease condensation into domains: 
results are only meaningful once this damped period has elapsed. 



3 Kinky vortons 

The kinky vorton model allows a charged condensate to live on a domain wall kink solution. 
The model has Lagrangian 

C = + d^a - ^ (0 2 - 4) 2 - ± (M 2 - r? a f - /^Vl 2 . (2) 

The first two terms are the kinetic terms of the real and complex scalar fields, (f>, a, respectively. 
Third is the symmetry breaking term of the real scalar - this allows domain walls to form. Next 
is a term which has a U(l) symmetry, where the complex scalar takes on a non-zero value in 
the vacuum. Finally, a quartic interaction term between the fields. This theory has a global 
Z2 x U{\) symmetry. The model parameters are chosen such that the Z2 symmetry is broken 
in the vacuum, and the U(l) symmetry retained. A consequence of the U{\) symmetry is that 
there exists a conserved symmetry current and charge (Noether's theorem): 

J M := ad^a - ad^a, fyj" = 0, — / d 3 x J = 0. 

One will immediately notice that there are a large number of model parameters, whose values 
must be chosen. A range of different choices have been discussed in ref.^, but in the present 
paper will only consider a single parameter set: = A CT = 2, /3 = 77^ = 1, r\ c = V3/2. 

A kinky vorton is a stable solution to the Euler-Lagrange equations of motion, which cor- 
responds to a circular kink solution with a current-carrying charged condensate living on the 
domain wall^. The radius of the loop is entirely determined by the charge and current that 
reside on the loop: one can imagine the loop being stabilized by the existence of the conserved 
charge/current (without the condensate, the loop collapses). 



3.1 Formation from random initial conditions 

The loops just described are entirely idealized: single loops, carefully setup in numerical simu- 
lations. So, we ask the question: can these charged loops form from random initial conditions, 
and can they aid network stabilization? 

We initially setup a homogeneous charged background, and random vacuum occupation. 
We specify an initial charge density /5q(0), and let the equations of motion evolve. Figure [Q 
has images of the 0-field, over time, for various initial charge densities /Oq(0). It is clear from 
inspecting these images that as the initial charge is increased, the rate at which the network 



"dilutes" slows down. This is quantified by inspecting a plot of the evolution of the total number 
of domain walls, for various initial charge densities (not included in the present paper, but can 
be found another publication^). We have also isolated some of the loops, and analyzed their 
properties (winding number, charge and radius). We have shown that they have properties very 
close to kinky vortons. 




Figure 1: Left Images of the evolution of the i^-field at t = 80, 160, 320, 640, 1280 (left to right) for initial charge 
densities Pq(0) = 0,0.01,0.09,0.25 (bottom to top). Right Image of 5R(ct), of the Pq(0) = 0.09 simulation at 
t = 640. In both sets of images, red/blue denotes positive/negative values (i.e. the two vacua, and a winding, 

respectively). 



4 Charge-coupled cubic anisotropy model 

The charged-coupled cubic anisotropy model is constructed by interacting a real vector field 
(f){x^) G R with a complex scalar field a, so that the model has a broken O(N) x U(l) global 
symmetry. It is only the O(N) symmetry that is broken in the vacuum - so that domain walls 
form, and there exists a conserved charge. The potential is 

\ N \ 

V = f (H 2 - 4? + £ £ 4 + ^ (M 2 - rff + m VI 2 - (3) 

1=1 

The first two terms create discrete vacua in the vector field, then the introduction of the con- 
densate field, and finally the term which interacts the condensate and vector fields. As there 
are more than 2 vacua, junctions and multiple tensions of domain walls can form (a junction 
is where more than two vacua cycle around a point in space: if four vacua cycle, an X-type 
junction is formed, and if there are three then a y-type junction). A known result is that the 
number of domain walls in the cubic anisotropy model drops n oc t : junctions do nothing to 
prevent the collapse of a network^. 

We have performed a set of simulations with this model, taking 4>{x^) G M. 2 . As in the kinky 
vorton model, we start off with homogeneous initial charge and random initial domain occupation 
(this time there are 4 vacua/colours to choose from). The parameter choice in this model becomes 
less trivial, but an interesting parameter worthy of note is the interaction constant /3, where we 
take 4/3 2 = A^Ao- + 4eA CT , chosen to avoid the phase-separation phenomenon^. 

In Figure [2] we display images of the 0-field, all of which are at the same time-step, but with 
different initial charge densities. Again, the implication is clear: more charge freezes the pattern 
in more. Plots of K(<r), charge and current densities look qualitatively identical to their kinky 
vorton counterparts, having become associated exclusively with domain walls. 




Figure 2: Images of the 0-field at t = 320, for p Q (0) = 0.09, 0.16, 0.2, 0.25, 0.36. 



5 Conclusions and summary 

Our first result is that networks with higher initial charge density freezes wall networks in 
"more" . Secondly, charge and current condense on domain walls, and almost exclusively live on 
them. Finally, some of the loops that form have the specific properties of kinky vortons. 

The elastic dark energy model was proposed to use domain walls as a model, requiring the 
wall network to freeze in. The conserved charge of the condensate field discussed here provides a 
resisting force against the tension of a domain wall, enabling a network to become pseudo-stable. 
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